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Abstract 

In this paper characterizations of graphs satisfying heat kernel estimates for a wide 
class of space-time scaling functions are given. The equivalence of the two-sided heat 
kernel estimate and the parabolic Harnack inequality is also shown via the equivalence 
of the upper (lower) heat kernel estimate to the parabolic mean value (and super mean 
value) inequality. 
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1 Introduction 

The heat propagation through a medium is determined by its heat capacity and conductance 
of the media. For the Euchdian space this observation goes back to Einstein. The heat 
diffusion has been a subject of interest in the discrete and continuous case for several decades 
and the fundamental results go back to the classical works of Aronson [1] , Davies [9] , Fabes, 
Stroock [11], Grigor'yan, [12], Moser [25], [26], Li,Yau [24], Saloff-Coste [29], Varadhan [39]. 
All these works are confined to homogeneous spaces. The diffusion in these spaces is typically 
located within the distance at time t, from the starting point. In other words, the time- 
space scaling is (time)^^^ or the space-time scaling is (distance)'^ . The inhomogeneous case 
attracts more and more attention of physicist and mathematicians since the 80-s. Geometric 
and algebraic conditions are relaxed and fractals enriched the topics. ( For recent results see 
[8], [13], [28], [20], [7], [4] and [5].) 

Dclmottc has shown in [10] (in the spirit of the results for manifolds by Saloff-Coste [29] 
and Grigor'yan [12]) for general graphs that the two-sided Gaussian heat kernel estimate 

c exp f-C^^) C exp (-c^^\ 

r\ <Pn{x,y)< " ^ (1.1) 

is equivalent to the parabolic Harnack inequality (with R = scaling, see all the formal 
definitions below). 

In the last two decades several works have been devoted to fractals and fractal fike 
graphs. One of the particular features of these structures is that the walk (or process) 
admits the space-time scaling function R^ with an exponent (3 > 2. For the continuous case 
the equivalence of the two-sided heat kernel estimate and the parabolic Harnack inequality 
with F {x, R) = R^ scaling has been shown in Hebisch, Saloff-Coste, [20] (see also Barlow, 
Bass [4] and Barlow, Bass and Kumagai [5]). In the graph case, the equivalence to the 
two-sided sub-Gaussian estimate (for (3 > 1) 



cexp 



-C (^) j Cexp [-C (^) 

< Pn {x, y) < 



V 



;i.2) 



and other conditions was shown in [15]. For a wider set of space-time scaling functions 
F [x, R) = F (R) the corresponding results were obtained in [36]. 

Barlow, Coulhon, Grigor'yan investigated in [6] the long time behavior of the heat 
kernel on manifolds using volume growth conditions. Here a more detailed picture will be 
provided covering on- and off-diagonal estimates under volume growth and potential theoretic 
conditions. 

Among others Hino, Ramirez [19], Norris [27] and Sturm [32] (see also references there) 
studied the heat diffusion in Dirichlet spaces. Their approach uses the intrinsic metric which 
recovers the classical Gaussian heat kernel estimate. For us the metric is a priori given and 
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the space-time scaling might be different from the classical R^. (For more comment about 
the difference of the two approaches with respect of fractals see comments in Section 3.2 of 
[19].) 

The present paper is partly motivated by the works Li and Wang [23] and Sung [30]. 
We prove in the context of weighted graphs that for a wide set of scaling functions the heat 
kernel upper estimate is equivalent to the parabolic mean value inequality (among others 
this is shown in [23] for the R'^ scaling). We also show ( inspired by [30] confined to the 
R"^ scaling) that some lower estimates are equivalent to the super mean value inequality. As 
a consequence, we prove that the conjunction of the parabolic mean value and super mean 
value inequality is equivalent to the two-sided heat kernel estimate and to the parabolic 
Harnack inequality, as well. 

Recent studies successfully transfer results obtained in continuous setting to the discrete 
graph case and vice versa (c.f. [5], [15], [16] ). For instance in [5] the proof of the equivalence 
of the parabolic Harnack inequality and two-sided heat kernel estimate (for the R^ scaling) 
is given for measure metric Dirichlet spaces via to the graph case where the equivalence is 
known (c.f. [4]). We treat the graph case while we believe that all the arguments and results 
can be transferred and are valid for measure metric spaces equipped with a strongly local, 
regular symmetric Dirichlet form and with the corresponding diffusion process. 

The aim of the present paper is to relax, as much as possible, the conditions imposed on 
the space-time scaling function. We will consider graphs for which the space-time scaling 
function F {x, R) is not uniform in the center x. One may feel that such a generalization 
is formal. A very simple example shows the opposite (see [36]), for the constructed graph 
neither the volume V {x, R) nor the space-time scaling function F [x, R) is uniform in x e F 
but heat kernel estimates hold. 

Fractafolds, defined by Strichartz [31], are the continuous counterparts of such struc- 
tures and as it is mentioned above we expect that the presented results are transferable to 
continuous spaces and to fractafolds. 

Before we can state our results we need some definitions. 



1.1 Basic definitions 

Let us consider a countable infinite connected graph F. A weight function /i^ y — ^y^^ > is 
given on the edges x ~ y. This weight induces a measure ^{x) 

y^x yeA 

on the vertex set A C F and defines a reversible Markov chain e F, i.e. a random walk 
on the weighted graph (F, //) with transition probabilities 

fJ,{X) 

Pn{x, y) = P(X„ = y\XQ = x) and the corresponding kernel, 
Pn{x,y) = -^Pn{x,y). 
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Let us use the notation p„ = Pn+i + Pn- The graph is equipped with the usual (shortest 
path length) graph distance d{x, y) and open metric balls are defined for x G F, i? > as 
B{x, R) = {y & r : d{x, y) < R}. The /x— measure of balls is denoted by 

V(x,R) ^ li(B(x,R)). (1.3) 

For a set A C F the killed random walk is defined by the transition operator restricted to 
Co (A) (to the set of functions with support in A) and the corresponding transition probability 
and kernel is denoted by {x, y) and (x, y). 

The (heat) kernel p„ (x, y) is the fundamental solution of the discrete heat equation on 

dnU = Am, (1.4) 

where dnU = m„+i — Un is the discrete differential operator with respect of the time and 
A = P — / is the Laplace operator on F. 

Definition 1.1 Throughout the paper we will assume that condition (Pq) holds, that is, there 
is a universal po > such that for all x,y & r,x y 

^>Po. (1.5) 

Definition 1.2 The weighted graph has the volume doubling (VD) property if there is a 
constant Dy > such that for all x & V and R> 

V{x,2R) < DvV{x,R). (1.6) 

Notation 1.1 For convenience we introduce a short notation for the volume of the annulus: 
v{x, r, R) = V{x, R) - V{x, r) for R > r > 0,x e F. 

Definition 1.3 Now let us consider the exit time 

TBi:c,R)^T^iHk-Xk^B{x,R)} 

from the ball B{x, R) and its mean value 

E,{x,R)^E{TBi,,R)\Xo^z) 

and let us use the notation 

E{x,R) = E^{x,R). 

Definition 1.4 We will say that the weighted graph (F,/x) satisfies the time comparison 
principle (TC) if there is a constant Ct > 1 such that for all x and R > 0,y G B {x, R) 
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Definition 1.5 We will say that the weighted graph (F, /i) satisfies the weak time comparison 
principle (wTC) if there is a constant C > 1 such that for all x & V and R > 0,y & B {x, R) 

#4 < C. (1.8) 

Notation 1.2 For a set ^4 C F denote the closure by 

A = {y eT : there is an x E A such that x y} , 
the boundary dA — A\A and A'^ — r\A. 

Definition 1.6 A function h is harmonic on a set A (ZT if it is defined on A and 

Ph{x) = J2Pi^,y)hiy)=h{x) 
y 

for all X E A. 

Definition 1.7 The weighted graph (F,//) satisfies the elliptic Harnack inequality (H) if 
there is a C > such that for all x eT and R> and for all u > harmonic functions on 
B{x, 2R) the following inequality holds 

max u < C mill u. (1-9) 

B{x,R) B{x,R) 

One can check easily that for any fixed Rq for all R < Rq the Harnack inequality follows 
from (po)- 

Definition 1.8 We define Wq to be the set of functions which are candidates to be a space- 
time scaling function. In particularF e Wq i/ F : F x N — > R and 
1. there are /3 > 1, > 0, c^, > such that for all R > r > 0,x E r,y e B {x, R) 



ci^Y^^^cJ^] , (1.10) 



\f{x^ ^^^^ 
r) ~ F{y,r) 

2. there is a c> such that for all x E T , R > 

F{x,R)>cR'^, (1.11) 

3. 

F{x,R + l)> F{x,R) + l (1.12) 

for all ReN. 

Finally F E Wi if F E Wq and (3' > I holds as well. 

Remark 1.1 We have by (1.12) that the function F{x,R) is strictly increasing from N to 
M in the second variable consequently it has the generalized inverse f : N ^ N : 

f {x, n) ^mm{REN: F {x, R) > n} . 

In the whole sequel f {x, n) is reserved for this inverse. 



5 



The function sets Wi C Wo will play a particular role in the whole sequel. 

Sometimes we will refer to the upper and lower estimate in (1.10) for x = y as the 
doubling and the anti-doubling property and in general, jointly we refer to them as doubling 
or regularity properties. 

Definition 1.9 We say that PH (F), the parabolic Harnack inequality holds for a function 
F if for the weighted graph (F, fi) there is a constant C > such that for any x ^ T, R, k > 
and any solution u > of the heat equation (1.4) on V = [k, k + F{x,R)] x B{x,2R) the 
following is true. On the smaller cylinders defined by 

1 1 

V- = [k + -F{x,R),k + -F{x,R)]xB{x,R) 

3 

andV+ - [k + -F{x,R),k + F{x,R)) X B{x,R) 

and taking (n_,x_) G , {n^,x^) G V^, 

d{x^, x+) < n+ — n 

the inequality 

Un_ix-) < CUn+{x+) 

holds, where we use the Un — Un + Un+i ■ 

Definition 1.10 Let us define the "volume" of a space-time cylinder D = [n, m] x (x, R) 

where m > n, R> by 

u (D) = [m-n]V {x, R) . 

Definition 1.11 We say that PMV5 (F) (the strong form of) the parabolic m,ean value 
inequality with respect to a function F holds on (F, fi) if for fixed constants < Ci < C2 < 
C3 < C4 < C5, < (5 < 1 there is a C > 1 such that for arbitrary a; G F and R > 0, using 
the notations F ^ F {x, R) , B = B {x, R) , V = [0, c^F] x B, V' = [c^F, cgF] x B {x, SR) , 
V'^ — [csF, C4F] X B {x, SR) for any non-negative Dirichlet sub-solution of the heat equation 

on V, the inequality 



holds. 

Definition 1.12 We will use PMV (F) if PMV5 (F) holds for 6 = 1. 

Definition 1.13 We say that (the strong form of) the parabolic super mean value inequality 
PSMV (F) holds on (F, /i) with respect to a function F if there is anO < e < 1 such that for 
any constants < ci < C2 < C3 < C/i < C5, C4— ci < e, there are 6, c> such that for arbitrary 
X G F and R > 0, using the notations F = F {x, R) , B = B {x, R) ,V = [0, C5F] x B, 
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(1.13) 



C 

maxM < ^ My)M (1-14) 



V'^ — [C3F, C4F] X B{x,6R), T> — [ciF, C2F] x B{x,SR) for any non-negative Dirichlet 
super-solution of the heat equation 



on T>, the inequality 

c 

mm 



holds. 

Definition 1.14 We introduce for A cT 



k=0 

the local Green function, which is the Green function of the killed walk and the corresponding 
Green's kernel as 

g''{y,z) = ^G^{y,z). 

Definition 1.15 The Green kernel may satisfy the following properties. There are c,C > 
and a function F such that for all x e T, R > 0, A = B {x, R) \B {x, R/2) ,B = B{x, 2R) 

max^^ (x, y) < c fj^'^^j , (1.16) 
yeA^ ^ - V {x, 2Ry ^ ' 

r,, , F(x,2R) , , 

mmg^(x,y) > c ,/, 'J 1.17 
yeA^ ^ ' - V {x, 2R) ^ ' 

If both inequalities hold this fact will he denoted by g (F) . 
1.2 Statement of the results 

The main results of the paper are the following. 

Theorem 1.1 // a weighted graph (r,/^) satisfies (po) o,nd (VD), then the following state- 
ments are equivalent. 

1. There is an F & Wq such that g (F) is satisfied, 

2. (wTC) and (H) hold, 

3. there is an F E Wq such that the upper estimate UE (F) holds: there are C, (3 > l,c> 
such that for all x,y eT, n > 



Pn [x, y) < —, — 77 TT exp 

V{xJ{x,n)) 



— c 



F{x,d{x,y))\^-^ 



n 



(1.18) 
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and furthermore the particular lower estimate PLE (F) holds: there are < c, 5 , e < 1 
such that for all x & F, R > 0, B — B {x, R) , d {x, y) < n A Sf {x,n) , n < sF {x, R) , 



Pn y) > 



;i.i9) 



V{x,f{x,n)y 

where f {x, n) is the inverse of F (x, R) in the second variable, 
4. there is an F e Wq such that PMV (F) and PSMV (F) hold. 

Further equivalent conditions will be given in Section 4. 

Remark 1.2 The off-diagonal lower estimate LE (F) which states that there are C,(3' > 
1, c > such that for all x,y E T, n > d {x, y) 



Pn {x, y) > 



V(x,f(x,n)) 



exp 



-C 



F{x,d{x,y))\W-i 



n 



;i.2o) 



can be obtained from (VD) and PLE{F) if j3' > 1 in (1.10) using Aronson's classical 
chaining argument. This indicates the possibility to obtain two-sided heat kernel estimate 
and necessary and sufficient conditions for it. 

Theorem 1.2 // a weighted graph (F, //) satisfies {po), then the following statements are 
equivalent: 

1. iVD) holds and there is an F E W\ such that g{F) is satisfied, 

2. there is an F E Wi such that the two-sided heat kernel estimate hold: there are C, (5 > 

> 1, c > such that for all x,y E T, n > d {x, y) 



exp 



-C 



F{x,d) \ W 
n 



exp 



< Pn {x, y) < C- 



— c 



F{x,d) \ l^' 



V{x,f{x,n)) 



(1.21) 



V{xJ{x,n)) 
where we write d — d{x,y), 

3. there is an F e Wi such that PMV (F) and PSMV (F) hold. 

4. there is an F E Wi such that PH (F) holds. 

Let us remark that the observation that the pair of the upper and particular lower 
estimate is equivalent to the /?— parabolic Harnack inequality goes back to [20]. 

Our results are presented in the discrete case, and with this limitation (which we consider 
not essential) are generalization of several works devoted to heat kernel estimates and the 
parabohc Harnack inequahty for scaling function F (x, R) = R^, R^ or F (i?), among others 
[29],[12],[8],[14],[15],[16],[30]. 

The following elements of the present paper are new: 
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1. the wide sets Wq and Wi of space-time scaling functions, 



2. condition g (F) with respect to F G Wi, i = 0,1, 

3. the parabohc inequahties with respect to F e Wi, i = 0, 1, 

4. the proof of the equivalence of the conjunction of the parabolic mean- and super mean 
value inequality to the parabolic Harnack inequality, 

5. the role of the strong anti-doubling property, > 1 is explained. 

The condition g (F) is the generalization of the corresponding conditions (G) in [14] and 



In [30] partial equivalence (for F (x, K) = B?) was shown for the parabolic super mean 
value inequality and the particular heat kernel lower estimate. Here we prove full equivalence 
for all F G Wq for a slightly modified version of the parabolic super mean value inequality. 
This modification allows us not only to show the full equivalence, but also appropriate for 
provcing that the conjunction of the parabolic mean and super mean value inequality is 
equivalent to the parabolic Harnack inequality. We have not found such a result in the 
hterature even for the classical case F [x, R) ~ R"^. 

In most earlier works the space- homogeneous case E {x, i?) ~ F (R) is considered. In 
these situations f3' > 1 follows from the homogeneity (and from the other conditions needed 
for the heat kernel estimates). 

The structure of the paper is the following. Section 2 contains the basic definitions. 
Section 3 recalls the results regarding the mean exit time. Section 4 contains the proof 
of Theorem 1.1, first the Subsection 4.1 summarizes a result about the heat kernel upper 
estimate, and Subsections 4.2 and 4.3 the lower estimate. The proof of Theorem 1.2 which 
contains the parabolic Harnack inequality is given in Section 5. The paper is closed with a 
short remark on the homogeneous case. 



(G^) in [15]. 



2 Definitions and preliminaries 



2.1 The volume 



Definition 2.1 We will use the inner product with respect to n : 



(/' 9)f, = ^f {x) g {x) // {x) . 



Remcirk 2.1 One can show that (VD) is equivalent to 




where a — log2 Dy and d{x, y) < R. 
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Remark 2.2 It is easy to show (c.f. [8]) that the volume doubling property implies the 
anti-doubling property: there is an Ay > 1 such that for all x & r,R> 

2V{x,R) <V{x,AvR) (2.1) 

which is equivalent with the existence of c,a'>0 such that for allx&r,R>r>0 

V(x,R) fR 
> c ' 



V{x,r) 

Notation 2.1 For two real series a^,b^,^ E S we shall use the notation ~ b^ if there is 
a C > 1 such that for all ^ & S 

Remark 2.3 Another direct consequence of (po) o,nd (VD) is that 

V (x, R, 2R) = V{x, 2R) - V(x, R) ~ V{x, R) (2.2) 

2.2 Laplacian 

Definition 2.2 The random walk on the weighted graph is a reversible Markov chain and 
the Markov operator P is naturally defined by 

pf{x) = Y,P{^.y)f{y)- 

Definition 2.3 The Laplace operator on the weighted graph (F, //) is defined simply as 

A = P-I. 



Definition 2.4 The Laplace operator with Dirichlet boundary conditions on a finite set A C 
r is defined as 

Afi^^_\ A/ {x) if xeA 
if x^A - 



A^/ (x) 



The smallest eigenvalue of —A"^ is denoted in general by X{A) and for A — B{x, R) it is 
denoted by X — X{x, R) — X{B{x, R)). 

Definition 2.5 The energy or Dirichlet form E {f,f ) associated to (r,/^) is defined as 

8 if, f) = - (A/, f), = \Y. fi.,y if {x) - f {y)f . 

x,yer 

Using this notation the smallest eigenvalue of — A"^ can be defined by 

A (A) = inf I ^jl^ : / e Co (A) , / ^ ^ (2.3) 

as well. 
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2.3 The resistance 

Definition 2.6 For any two disjoint sets, A,B(ZV, the resistance, p{A,B), is defined as 

p{A, B) = (inf {£ if, /) : = 1, /l^ = O})"^ 

and we introduce 

p{x,r,R) ^ p{B{x,r),r\B{x,R)) 
for the resistance of the annulus around x e F, with R > r > 0. 

Definition 2.7 We say that the product of the resistance and volume of the annulus is 
uniform in the space if 

p{x, R, 2R)v{x, R, 2R) ~ p{y, R, 2R)v{y, R, 2R). (2.4) 

We will refer to this property shortly by (pv) . 

Lemma 2.4 For all weighted graphs, xer,R>r>0 

p{x, r, R)v{x, r, R) > {R - r f, (2.5) 

Proof. For the proof see [35]. ■ 

Definition 2.8 The resistance lower estimate RLE (F) holds for a function F if there is a 
c > such that for all x & r,R> 

Definition 2.9 The anti-doubling property (aDpv) is satisfied for pv if there are c, (3' > 
such that for allx&r,R>r>0 

p {x, R, 2R) V (x, i?, 2R) ^ ^ fR\ 



p (x, r, 2r) V (x, r, 2r) \ ^ / 

2.4 The mean exit time 

Let us introduce the exit time Ta from a set ^4 C F. 

Definition 2.10 The exit time from a set A is defined as 

Ta = min{A; : Xk G F\A}, 

its expected value is denoted by 

E,{A) = E(7:4|^o = x) 

and furthermore leus us write 

E ix, R) = max Ey (B (x,R)) 
11 



In this section we introduce some properties of the mean exit time which will play crucial 
role in the whole sequel. First of all it is immediate that 



and we for e N 



E {x, 1) > 1 



E{x,R + l) >E{x,R) + 1. 



(2.8) 



Remcirk 2.5 We have by (2.8) that the function E {x, R) is strictly increasing from N toW 
in the second variable consequently it has the generalized inverse e : N — > N; 

e {x, n) — min {R E N : E {x, R) > n} . 

Remark 2.6 It is easy to see that (TC) is equivalent to the existence of constants C, /3 > 1 
for which 



E{x,R) 



< C 



(2.9) 



for ally e B{x,R) ,R> S > 0. 



Definition 2.11 The local sub-Gaussian upper exponent, with respect to a function F {x, R) 
is k — kx{n, R) > 1, it is defined as the maximal integer for which 



n 



— < q min F(y, 

k y&B{x,R) 



R 



) 



(2.10) 



or k = 1 by definition if there is no appropriate k. Here q > is a small fixed constant (c.f. 
[34],q < min {1/16, CfPo/Cf};. 



n 



Definition 2.12 Let n > l^ — Ixin, R) > 1 be the minimal integer for which 

) (2.11) 

or I = n by definition if there is no appropriate I. The constant C will be specified later. 



J > CF{x, 



R 
J 



Definition 2.13 The local sub-Gaussian lower exponent I {n, R, A) with respect to a function 
F {x, R) for A (ZT is the maximal integer I for which 



n 



— > C max F(z, 

I z€A ^ 



R 
1 



(2.12) 



Definition 2.14 The global sub-Gaussian exponent m = m{n, R) is defined as the maximal 
integer for which 

— < qminEiy, — ) (2.13) 
or m — 1 by definition if there is no appropriate m. 



12 



Definition 2.15 The mean exit time is uniform in the space if there is a function F such 
that 

E{x,R)c^F{R). (2.14) 
This property will be referred to by {E) . 

The definition m (n, R) is prepared for the particular case when E {y, R) c::^ E [x, R) ~ 
F (R) i.e. E is basically independent oi x,y e F. 

Definition 2.16 We define V\ to be the a set of functions such that F E Vq if F : FxN ^ M 
and there are (3' > 1,cf > such that for all R > r > 0,x E T,y E B {x, R) 

^\r) - F{y,r) ^ > 



2.5 Mean value inequalities 

Definition 2.17 The elliptic mean value inequality (MV) holds if there is a C > such 
that for all X E T , R > and for all u >0 harmonic functions on B = B (x, R) 

^ ' ' y&B 

Remark 2.7 Let us recognize that in the defimtion of PMV S is a "free" parameter, while 
in the definition of PSMV {F) it depends on e and Ci,i = 1..5. Let us also observe that Ci-s 
are subject of the restriction C4 — ci < £ . 

Remark 2.8 One should note that the definition of the parabolic mean value inequality is 
slightly different from that is given in [36]. There it is stated for Dirichlet solutions, here 
we have it for arbitrary Dirichlet sub- solutions. It is easy to see that the extended definitions 
fits into Theorem 4- 2- On one hand solutions are sub-solutions, on the other hand, the proof 
of the implication UE =^ PMV of Theorem 4-2 follows word by word for sub-solutions. 

Remark 2.9 The condition (1-13) in the definition of the parabolic Harnack inequality is 
needed in order to have a path (with nonzero probability ) of length no more than n+ — n_ 
between X- and x^. One can eliminate this restriction if the parabolic Harnack inequality 

is considered only for large enough R-s. The condition n+ — > d{x^,x^) is satisfied if 
(c-i — C2) F [x , R) > cR^ > 4R which holds if R > Rq. Such an Rq depends only on the 
constants (c.f [35]). In order to avoid lengthy technical discussion we may assume R > Rq 
in all these situations. If other is not stated, the corresponding inequalities for R < Rq 
follow from (po) • 
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3 Properties of the mean exit time 



In this section we recall some results from [35] which describe the behavior of the mean exit 
time. If this is not mentioned otherwise, the statements and proofs can be found in [35]. 
The first one is the Einstein relation. 

Theorem 3.1 // (po) , (VD), (H) and one of the conditions (wTC) , (aDpv) , RLE {E) or 
pv e Wq hold, then (ER), the Einstein relation 

E{x, 2R) ~ p{x, R, 2R)v{x, R, 2R) (3.1) 

holds, furthermore 

A"^ {x, R)^E {x, R)c^E (x, R) 

and 

E (x, R) e Wq. 

Theorem 3.2 // (po) , (YD), (TC) hold, then the Einstein relation 

E{x, 2R) ~ p{x, R, 2R)v{x, R, 2R) (3.2) 

holds, furthermore 

(x, R)c^E {x, R)c:^E (x, R) e Wo. 

The properties of the inverse function e and properties of E are linked as the following 
evident lemma states. 

Lemma 3.1 The following statements are equivalent 

1. There are C,c> 0, ^ > > such that for all x e F, R > r > 0, y e B {x, R) 



2. There are C,c>0,(3>(3'>0 such that for all x G T,n > m > 0, y ^ B {x,e {x, n)) 

V/3 e {x, n) / n \ V/?' 
e (y, m) 



c(^y'<^^<c(^y\ (3.4) 



Lemma 3.2 If F E Wq, then for kx{n, R) defined in (2.10) 

k^{n, R) + l>c (^^^^^^ " ' (3.5) 

for all X &r,R,n > for fixed c > 0, (3 > 1. In addition, if F & Wi is assumed, then 

k {n,R)-l<C ^ ■ (3.6) 
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Proof. The statement follows from the regularity properties of F easily, (3 > 1 is ensured 
by /? > 2 (see (1-11)) and > 1 by the assumption. ■ 

Lemma 3.3 // (E) is satisfied and E e Wq then 

kx (n, R) ~ Ix (n, R) m (n, R) . 

Since this fact is not used in the proof of the main results, the elementary proof is omitted 
(for some hints see [34]). 

Remcirk 3.4 Let us mention here that under (po) , (VD) and (H) the uniformity of the 
mean exit time in the space: 

E {x, i?) ^ F (R) 

ensures that E satisfies the left hand side of (1.10) with a f3' > 1 (c.f. [35]). This explains 
that in the "classical" cases, when (E) holds one should not assume j3' > 1, it follows from 
the conditions (see also [36] Theorem 4.13). 

Lemma 3.5 For (r, //) for all x er,R > 

min E {z, R/2) < p (x, i?, 2R) v (x, R, 2R) . (3.7) 

Corollary 3.6 Under (po) and (VD) 

3F eWo-.g (F) ^ (H) + (ER) . 

Proof. The implication was shown in [35]. We also know that (po) , i^D) , (H) and 
{ER) implies (TC) hence by Theorem 3.2 E G Wq. The reverse implication needs some 
additional arguments. We know again from [35] that g {F) implies {H). We show here the 
imphcation g (F) =^ (ER) under (po) , (VD) and (H). That needs some care. Let us assume 
that ri = 2\ r„_i < 2R < r^, Bi = B {x, r-) , , A = Bi\Bi^i, Vi^V {x, ri) . In [15] Section 
4.3 it is derived using (po) , {VD) and {H) that 

n— 1 

E (x, 2R)<CY, Vi+iP i^, n+i) . 

1=0 

Now we use a consequence of {H): 

p(x,rj,rj+i) < C max g'^'+^ {x,y) 

(see for instance [35] Section 4. or [3]) to obtain 

n— 1 

E{x,2R) < C^Vi+i max g^^+^ [x, y) 
^ 2/eAi+i 

1=0 

n— 1 n— 1 

< C 5^ F (x, r,+i) < CF (x, r„) 2"*^' 

1=0 1=0 

< CF{x,2R), 
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where (1-16) was used to get the second inequahty. 
On the other hand, from (1.17) one obtains 

^ ' « I. - M V , / ; zeB{x,R)\B{x,R/2) 

< Yl 9''{x,z)ii{z)<E{x,2R), 

z&B(x,R)\B{x,R/2) 

this means that 

cF {x, 2R) < E {x, 2R) 

consequently, F c:^ E, E & Wq and (TC) is satisfied. Finally by Theorem 3.2 the conditions 
(po) , (VD) and (TC) imply (EE), m 

Corollary 3.7 Assume that (F,//) satisfies (po) , (VD) and (H), then 

(wTC) ^ {aDpv) ^ {TC) ^ {EE) ^ ELE {E) ^ g {F) , (3.8) 



where F e Wq is a consequence in the direction =^ and assumption for 




Proof. Exept the last implications the statement was shown in [35] while the last one is 
just Corollary 3.6. ■ 

Remcirk 3.8 Let us remark here that as a side result it follows that ELE {E) or g {F) for 
F e Wq implies pv :^ F and E :^ F as well. 

4 Temporal regularity and heat kernel estimates 

In this section we prove Theorem 1.1 in an extended form. We have seen in Corollary 3.7 
that under the conditions (po) , {VD) and {H) 

{wTC) ^ {aDpv) ^ {TC) <^ {EE) <^ ELE {E) . (4.1) 

Let (*) denote any of the equivalent conditions. Using this convention we can state the 
extension of Theorem 1.1 as follows. 

Theorem 4.1 // a weighted graph (F, /x) satisfies (po) o,nd {VD), then the following state- 
ments are equivalent: 

1. there is an F & Wq such that g {F) is satisfied, 

2. {H) and {*) hold, 

3. there is an F & Wq such that U E {F) and PLE {F) are satisfied, 

4. there is an F e Wq such that PMV {F) and PSMV {F) are satisfied. 
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The proof of Theorem 4.1 contains two autonomous results. The first one states that the 
upper estimate is equivalent to the parabolic mean value inequality, the second one states 
that the particular lower estimate is equivalent to the parabolic super mean value inequality. 
The return route from 5 to 1 and 2 is based on the Einstein relation, on a potential theoretic 
results from [35] and a modification of the return route developed in [34]. The proof of 
2 =^ 3 generahzes methods of [14] and [15]. 

Let us emphasize the importance of the condition (aDpv) in (4.1). It is a condition on 
the volume and resistance, no assumption of stochastic nature is involved so the result is in 
the spirit of Einstein's observation on the heat propagation. These conditions in conjunction 
with (VD) and {H) provide the characterization of the heat kernel estimates in terms of 
volume and resistance properties. Of course the elliptic Harnack inequality is not easy 
to verify. Meanwhile we learn from g (F) that the main properties ensured by the elliptic 
Harnack inequality are that the cquipotential surfaces of the local green kernel g^(^'^^ are 
basically spherical and the potential growth is regular (c.f. [35]). 



4.1 The upper estimate 

This section provides the upper bound part of the implication 2 =^ 3 of Theorem 1.1. In 
details 

(VD) ] 

(TC) \ ^ DUE (E) (4.2) 
(H) J 

and under (VD) and {TC) 

DUE (E) UE (E) PMV (E) 

In particular the parabolic mean value inequality is shown to be equivalent to the upper 
estimate and the other conditions. This result has been proved in [36]: 

Theorem 4.2 For a weighted graph (F,//) if (po) , (VD), (TC) conditions hold, then the 
following statements are equivalent: 

1. The local diagonal upper estimate DUE [E) holds; there is a C > such that for all 
X er, n> 

C 

Pn{x,x)<— --, (4.3) 

V[x, e[x, n)) 

2. the upper estimate UE{E) holds: there are C, /3 > l,c > such that for all x,y &r, 
n>0 



Vn y) < — exp 

V [X, e [x, n)) 



—c 



E{x,d{x,y))\^ 



n 



3. the parabolic mean value inequality, PMV (E) holds, 
4- the mean value inequality, (MV) holds. 
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Corollary 4.1 //(F, //) satisfies (po), then 

=^ UE (E) . 



iVD) 
(wTC) 



Proof. We know from Theorem 3.2 that E e Wq. The statement follows from Corollary 
3.7 and Theorem 4.2 since the elliptic Harnack inequality, (H) implies the elliptic mean value 
inequahty {MV). m 

Remcirk 4.2 The equivalences in (??) for F e Wq instead of E follow from the same proofs 
given in [36] for E (see also Corollary 3.7 and Remark 3.8) 

Remcirk 4.3 Conequently the implication of the upper bound part in Theorem 4-1 (and 
Theorem 1.1) 2. ^ 3. is shown. 

Remcirk 4.4 The implication for an F e Wq 

DUE (F) ^ UE (F) ^ PMV (F) 

in particular DUE [F] =>■ PMV {F) can be shown repeating step by step the proof given for 
the particular space-time scaling function E. The full proof is spelled out in /?, Theorem 
8.6]. This gives the proof of the upper bound part of the implication of Theorem 4-13.^ 4. 

4.2 The near diagonal lower estimate 

In this section we give a lower estimate for the Dirichlct heat kernel and for the global one. 

Definition 4.1 The near diagonal lower estimate, NDLE (F) holds with respect to a func- 
tion F if there are c,S>0 such that for all x,y & T,n > 0,d {x, y) < Sf {x, n) An 

Pn {X, y) > TTJ—^i TT, (4.4) 

V [x,f [x,n)) 

where f is the existing inverse of F in the second variable, defined in Remark 1.1. 

Remark 4.5 R is clear that PLE (F) implies NDLE (F). R is also known (c.f.fl3]) that 
NDLE (F) and UE (F) implies PLE (F) if F ^ R^ , the same proof works F e W^o- 

Theorem 4.3 For weighted graphs 

(po) + (VD) + (TC) + (H) =^ PLE (E) 

The proof closely follows the steps of the corresponding proof given for the case E [x, R) ~ 
R^ in [14] therefore it is omitted. 

Remark 4.6 From the regularity ofV and F (and f) it is immediate that PLE (F) is equiv- 
alent with the slightly stronger form with 6' = 6/2 : for all x & T, R > 0, B = B [x, R) ,n < 
eF (x, K),y,z eB {x, 6'f {x, n)) 
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4.3 The parabolic super mean value inequality 

In this section the equivalence of the particular lower estimate and a kind of converse of 
the parabolic mean value inequality is shown. The partial equivalence for the classical 
{F {x, R) = R^) and continuous situation was shown in [30]. Here the generahzation to the 
present settings is provided. 

For technical reasons we use some specific constants, like £, 5 from PLE, cp, Cp from the 
definition of the set of scahng functions Wq (in (1.10)). 

In this section we show that for an F e Wq 

PLE (F) <^ PSMV (F) 
that is, the following theorem holds. 

Theorem 4.4 For the weighted graph (r,/i) assume (po) Q-nd {VD). Then for an F & Wq 
PLE (F) holds if and only if PSMV (F) holds as well. 

Proof. The proof follows the main steps of [30]. We know that PLE (F) implies the 
slightly stronger version (4.5) that is there are S,e > and c > 

pS(!/.^)+pS«(!/.-)>^7^^^^ (4.1) 

holds, provided that y,z E B {x,r) , where r = 6f{x,m)/2, and m < eF = eF{x,R). For 
the super-solution u we have that for all C3F < k < C4F, ciF < i < C2F 

Uk{y)> Yl Pk-i{v^^)l^{z)ui{z). 

zeB(x,R) 

In order to use (4.5), we choose 

S*^^{cs- 02)^5/2 (4.2) 

which ensures that y,z E B {x, r) if r = 6*R. From the condition C4 — ci < £ it follows that 
k — i < sF [x, R) is satisfied and PLE (F) can be applied: 

y&B{x,6*R) 

Now let us sum for ciF < i < C2F and divide by (c2 — Ci) F to obtain 

C2F 

- v(x R)F(x R) ^ ^ ^^iy)^^iz)■ 

Now we prove the reverse implication PSMV (F) =^ PLE (F) by applying PSMV 
twice. 
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1. Denote e, 5, Ci the constants in PSMV (F), q will be specified later (which determines 
5 as well), furthermore F = F {p, R) , B = B {p, R) , n = R/8, Fi = F {p, ri) , = 
B {p, Sri) , ^ = c'Fi, Ci < c' < C2 and D = B (j),S^y Let us define 



Un (y) 



J2zeDPn-m{y:^)l^{z) if 71 > m 

1 if n < m 

This is a solution on D x [0, oo] of 

and Un > 0. Prom the PSMV (F) it follows that 

C2F1 

'''^''^- v(xt,)F ^ E/^H^^H 

provided that, x e Di and C3F1 < k < C4F1. Prom the definition of ii„, (VD) and 
F e Wo it follows that 

C2F1 

Again from the definition of Un and Di C C3F1 < k < C4F1, x e Di we obtain 

J]pf_^(a;,^)/.(^)>c (4.3) 

or equivalently 

J]pf (a;,2;)/i(2;) >c (4.4) 

if X e (C3 - C2) Fi < i < (C4 - Ci) Fi. 

2. We will use the parabolic super mean value inequality in a new ball B2 for pf {x,y) 
with the same set of constants q, hence with the same 5 as well. Let r2 = -R/2, 
B2^ B {x, rs) , D2 = B {x, 5r2) , F2 = F (x, ra) . We apply PSMV in and obtain 
that for C3F2 < / < C4F2, y e -D2 

c '^^^^ 
V a;, dr2 F r2 

if in addition B2 C B {p, R). Let x & B {p, This ensures that B2 C B {p, R) and 
D2 D D and we obtain for y e S (x, ^) C F'2, (and B (x, ^ ^ (P' ^) as well) 
that 

C2F2 

^ y(x,^r2)F(a;,r2) ^ ^ " ^'"'^ 
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In order to use (4.4) we require 

C2F2 > (C4 - ci) Fi (4.6) 

and 

C1F2 < (c3-C2)Fi. (4.7) 
From the assumption F G Wq it follows that (4.6) is satisfied if 



C4 = C2 



(1 + cAf"') 



and (4.7) is satisfied if 



for any < g < 1. Finally 



(C3 - C2) B 

ci = q — 4 " 

Up 



< C4 - ci = C2 (1 + cj.4^') - gfe^^-M4-/3 < s 

and ci < C2 < C3 < C4 can be ensured with the appropriate choice of C2, C3 and 5. Using 
(4.6) and (4.7) and D2 D D the estimate in (4.5) can be continued as follows: 

pf {x, y) 

C2F2 

^ V{x,Sr2)Fix,r2) ^^^^'^y^'^^^'^- 

^ (c4-ci)Fi 

l-(x,5r2)F(x,r2) (^^t;')^^,^ 



Now we apply (4.4) to conclude to 



(c4-ci)Fi 



5^^'^^ (x v) > - V 



y(a:,5r2)F(a:,r2) ' 

c n 
> > 



V{x,R) - V{x,f{x,l)y 
where C3F2 < I < C4F2 and y e S (p, ^). Finally let 5 > 2i? 

pf ^^'^^ (x, I/) > pf {x, y) > — — ^- — — (4. 



under the same conditions. Now choosing e' — ^2 ^ and S' = ^ i^sCf)^^^ (4.8) 
implies that PLE{F) (in the stronger form: (4.5)) 



V{xJ{xJ)) 
^2-/^ and 5' - , 



foTdix,y)<5'fip,l),l<e'F{p,S) 
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4.4 Time comparison 

In this subsection we summarize the results which lead to the proof of 1 <^=^ 2 =^ 3 4 
in Theorem 4.1 and we prove the return route from 4 ^> 2 The equivalence of 1 and 2 is 
established by Theorem 3.2, Corollary 3.6 and 3.7, see also Remark 3.8. The implication 
2 =^ 3 is given by Theorem 4.2 and 4.3 and 3 =^ 4 is combination of Theorem 4.2 and 4.4, 
see also Remark 4.4. 

Now we prove 4 =^ 2, the return route of Theorem 4.1. Our task is to verify the 
implications in the diagram below under the assumption F e Wq and (po) , (VD). 



PMVi{F) \ PMVs*{F) 
PSMV (F) J ^ PSMV (F) 

DUE{F) 



(H) (4.10) 



PMV^iF) \ ^ { I pv^F\ [TC) 

P5My (F) / ^ ^^P^ W ^ {H) {H) ^^-^^^ 

The heat kernel estimates are established as we indicated above. Now we deal with proof 
of the elliptic Harnack inequality [H) and the time comparison principle {TC). 

Theorem 4.5 // V satisfies (po) , iVD) and there is an F e Wq for which PMV (F) and 
PSMV (F) are satisfied, then the elliptic Harnack inequality holds on (H). 

We need an intermediate step, the parabolic mean value inequality for smaller balls. We 
choose a particular set of constants q subject some restrictions coming from PSMV and 
needed for later use. 

Lemma 4.7 // (F, /j.) satisfies (po) , (VD) and PMVi (F) for an F e Wq, then for a given 

1 . 

e,S > 0,0 < S* < ^£^1 there are ci < ... < C4 such that PMV^* (F) holds for e and Ci — s. 

Proof. We would like to derive PMVg* (F) for q from PMVi (F) which holds for some 
other constants Oj. We will apply PMVi (F) on the ball B — B [x, 5R) and re-scale the time 
accordingly. We have PMVs* (F) on B {x, R) by 

max Ui (y) < max Ui (y) 

C3F{x,R)<i<C4F{x,R) a3F{x,SR)<i<a4F{x,SR) 
yeB yeB 



^ a2F{x,5R) 

^ ^ J]%(^)//(^) 

^ ' j=aiF(x,SR) yeB 
^ C2F{x,R) 



j=ciF{x,R) yeB 



if the inequalities ci < ... < C4, oi < . 


.. < 04, 








a^F 


{x, SR) 


> 


C4F 


{x,R) 


a^F' 


[x, 6R) 


< 


C3F 


{x,R) 


a2F 


[x, 6R) 


< 


C2F 


{x,R) 


aiF 


[x, 6R) 


> 


ciF 


{x,R) 
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1 „ 

are satisfied. We require in addition that C4 < e and ^ (C3 — 02)^ | > S*. One can see 
that the following choice satisfies these restrictions. Denote p — Cp {5*)^ ,q — cp (S*)^ . Let 



C4 




E, G4. — 


2g 
— C4 




C3 


< 


C4, 03 = 


1 

: - min 
2 




C2 


< 


ca, a2 = 


{PC2, 03} 


Cl 




1 . J 
- mm < 

2 1 




1 , ai = gci 



Let us observe that ci can be arbitrarily small since C4 < £ and if the sub-solution is not 
given from an m up to a^F (x, SR) it can be extended simply by Ui^m = Pf^^'^'^Um- ■ 

Proof of Theorem 4.5. Let us fix a set of constants ci < C2 < C3 < C4 = e as 
in Lemma 4.7 and apply P5'My (F) for them. Let us apply Lemma 4.7 for S* to receive 
PMVs* (F) on B — B (x, i?) . As a consequence ior D — B {x, S*R) , Uk (y) = h [y] we obtain 

max/i < C (4.12) 

Similarly PSMV (F) yields 

miuh>c^h{y). (4.13) 
yen 

The combination of (4.12) and (4.13) gives the elliptic Harnack inequality for the shrinking 
parameter 5* . Finally {H) can be shown using the standard chaining argument along a finite 
chain of balls. The finiteness of the number of balls follows from volume doubling via the 
bounded covering principle. ■ 

Theorem 4.6 // (^0)7(^-0) hold furthermore there is an F ^ Wq for which PMV (F) and 
PSMV {F) are satisfied, then E ^ F and {TC) is true. 

Proposition 4.8 Assume (po) and {V D) hold. If PLE {F) holds for F ^ Wq hold, then 
then there is a c> Q such that for all R> 0,x 

E {x, R) > cF {x, R) . 

Proof. It follows from PLE (F) that there are c, C, 1 > 5 > 5' > 0, 1 > e > > 
such that for all x e T, R > 1, A ^ B{x, 2R) and n : e'F {x, R) < n < eF {x, R) , r ^ 5'R, 
y e B ^ B {x,r) 

Pnix,y) = P:^{x,y) + P.^^,{x,y) > ^^^y 

It follows for F ^eF {x, R),F' ^ e'F {x, R) that 

00 00 

E{x,2R) = Y, p^{x,y)>Y^Y.2^k{^^y^ 

k=0 yeB{x,2R) k=0 yeB 

k=F'y£B ^ ' ^ 
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Proposition 4.9 // (po) , (VD) hold and DUE (F) holds for an F e Wq, then there is a 
C > such that for all R > 0,x eV 

p (x, 2R) V (x, 2R) < CF (x, 2R) . 

The first step towards the upper estimate of pv is to show an upper estimate for A~^. 

Proposition 4.10 // (po), {VD) ,DUE{F) hold and F e Wq, then there is a c > such 
that for all R > 0,x E r 

\{x,R) >cF-\x,R). (4.14) 
Proof. Assume that Ci > 1, 2n = \F{x, CiR)~\ , y,z e B = B(x. R). One can use 

-P2n(y, z)^^Pn (y, W) {w, z) < ^J Pan {V , V) P2n {z, z) 
w 

and DUE (F) to get 

p{z) 



P2n{y,z)<C 



{V{y,fiy,2n))V{z,f{z,2n))f'^ 



(for the details see [14]). Prom (VD) and F e Wq it follows ior w — y or z, d{x,w) < R < 
CiR = f{x, 2n) 

V{x, C,R) ^ 
V{w,CiR) - ' 

which results using (po) that for all n 

p{z) 



Pn{y,z)<C 



V{x,f{x,n))' 



If (f) is the left eigenvector (measure) belonging to the smallest eigenvalue A of — normal- 
ized to (01) — 1, then 

(1-A)- = 0P-1= Y: '^(-«(-'^)< E — ^v(l Hz 2.)) 

< C max „ . „ . = C max ' ^ ' ^ ' 



zeB{x,R) \f (z, 2n) J zeB{x,R) \Ci f (z, 2n) 

if Ci = 2CV"C/. Using the inequality and 1-^ > i log i for ^ e [|, 1], where ^ = 1-X{x, R), 
one has 

X{x,R) > > cF(x,CiP)-^ > cF{x,R)-\ 

■ 

Proof of Proposition 4.9. Let us recall from [35] that 

A {x, 2R) p {x, R, 2R) V {x, R)<1 
in general, applying {V D) and (4.14) immediately yields the statement. ■ 
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Proposition 4.11 // (po) , (VD) hold and PLE (F) for an F e Wq, then there is a c> 
such that for all R > 0,x & F 

p {x, R, 2R) V (x, R, 2R) > cF {x, 2R) 

Proof. The inequality (3.7) states tliat 

p {x, R, 2R) V {x, R, 2R) > min E (z, R/2) . (4.15) 

Prom Proposition 4.8 we know that 

min E {z, R/2) >c min F {z, R/2) 

zedB(^x,^R) zedB(^x,^R) 

and from F G Wq it follows that 

p {x, R, 2R) V {x, R, 2R) > min F {z, R/2) > cF {x, 2R) . 

z&dB{x,lR) 

■ 

Proof of Theorem 4.6. From Proposition 4.9 we have that pv < CF which together 
with Proposition 4.11 yields that 

p {x, R, 2R) V {x, R, 2R) ~ F {x, 2R) . 

Since F e Wo we have that pv e Wq and (aDpv) as well. Prom the conditions we have (H) 
thanks to Theorem 4.5 and by Theorem 3.2 the Einstein relation follows: 

E [x, 2R) ~ p (x, R, 2R) v [x, R, 2R) ~ F {x, 2R) . (4.16) 

Since F e Wq and ~ F it follows that E e Wq which includes (TC) and of course (wTC) 
as well and the proof of 4 =^ 2 of Theorem 4.1 is complete. ■ 

5 The parabolic Harnack inequality 

In this section we will prove the following extension of Theorem 1.2. 

Theorem 5.1 // a weighted graph (F, /x) satisfies (po), then the following statements are 
equivalent: 

1. iVD) hold and there is an F & Wi such that g{F) is satisfied, 

2. {VD) , (H) and {*) holds furthermore E e Vi, 

3. (VD) , (H) and pv e Vi, 

4. (VD) and UE (F) , PLE (F) , for anF eWi are satisfied, 

5. {VD) holds and there is an F & Wi such that PMV {F) and PSAIV {F) are true, 
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6. there is an F E Wi such that the two-sided heat kernel estimate hold: there are C,(3 > 
(3' > l,c> such that for all x,y & F, n > d {x, y) 



exp 
c- 



-C 



< Pn {x, y) <C 



(5.1) 



V(x,f(x,n)) -""^ V(x,f(x,n)) 
where d = d{x,y), 
7. there is an F e Wi such that PH (F) holds. 

The equivalence of the statements 1 — 5 are based on Theorem 4.1. What is left is to 
incorporate 6 and 7. In this section we show that the mean value inequahties for F e Wi 
are equivalent to the parabolic Harnack inequality and to the two-sided heat kernel estimate 
(5.1). We will show the following implications: 

PMV ] ( DUE ] ( (VD) 

PSMV \ ^ PH ^ I DLE \ =^ \ PMV , 
iVD) J 1 PSMV J I PSMV 



UE 
PLE 
(VD) 



UE 
LE 



Theorem 5.2 Assume (po). Let F e Wi, then the following equivalence holds: 

(VD) + PMV (F) + PSMV (F) <^ PH (F) . 

RemEirk 5.1 We give direct proof of the statement instead of the ready alternative from 
[10]. In Theorem 3.10 of [10] the decomposition method shows that for the parabolic Harnack 
inequality it is enough to show UE (F) and PLE (F) for the Dirichlet heat kernel on B (x, R) . 
Since we know that PMV <^==^ U E and PSMV <^==^ PLE the proof is similar hut works 
via the Dirichlet heat kernel estimates. Here we prefer the direct route. 

Proof. The proof consist of several smaller steps. 

1. First we show PH (F) for Dirichlet solutions for a particular set of constants. We 

1 

choose Ci < ... < C4 and = c> ('^3 ~ ^2)^ 5/2 as in Lemma 4.7. Denote = [C3F, C4F] x 
B (x, 5*R) and $" = [ciF, caFj^x B {x, 5*R). Using Lemma 4.7 we have for 8\ PMVg* (F) : 

C 

maxu< —^^^Ui{z)fi{z) (5.2) 

Let us choose cq > C5 > C4. The parabolic super mean value inequality PSMV (F) with 
V+ = [C5F, cqF] X B {x, S*R) , = states that 

min u > — ^ J2 (^) (^) ■ (5-3) 
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The combination of (5.2) and (5.3) results that 



maxM<Cminti (5.4) 
D- D+ 

which is the parabolic Harnack inequality for Dirichlet solutions for the constants C3 < C4 < 
C5 < C6, 5*, in other words D' = ^+ = V+ . 

2. Let us use the decomposition for an arbitrary solution -u; > on P = [0, F (x, R)] x 
B {x,R). The nonnegative linear decomposition results in a Dirichlet solution -u > on "D 
for which u = w on B [x, 6*K) and m < w in general (for the details of the decomposition 
method see [10] proof of Theorem 3.10). Now we use (5.4) 

max w — max u < C min u < C vainw. 

D- D- D+ D+ 

Which means that we have PH (F) for all solutions and for the given q— s and 6*. 

3. It is standard knowledge that if the (classical) parabohc Harnack inequahty holds for 
a set of constants q, S, then it is true for arbitrary set of constants as well (with an other 
C). This is the case a F e Wi. The key is that (3' > 1 ensures that the time dimension 
of the space-time cylinder shrinks faster than the space dimension and the usual chaining 
argument can be applied. 

4. The imphcation PH (F) =^ (VD) can be seen along the lines of the classical proof 
(c.f. [10]). First from PH (F) the diagonal upper and lower estimates are deduced without 
change of the proof 

C 

Pm (x, x) < — — — (5.5) 

V{xJ{x,m)) 



and 



The inequality for n < cm 



Pn {x, x) > — — — --. (5.6) 

V{x,f{x,n)) 



Pn {X, X) < Cpm {X, X) (5.7) 

can be obtained from PH (F) with the proper choice of the constants. Now let n = 

[F (x, R)\ ,m = \F (x, A^R)] , p > 1 and A > 2 is chosen to satisfy p > and A > 

Cp^Whs^ As a result from (5.5) , (5.6) and (5.7) one obtains (VD) : 



CF 



V {x, 2R) < V {x, APR) < CV {x, R) . 

5. The implication PH (F) =^ PSMV (F) is evident. As in step 4 we deduced 
PH (F) =^ DUE (F) and PMV (F) follows from Theorem 4.1. ■ 

Remcirk 5.2 The elliptic Harnack inequality is a direct consequence of the F-parabolic one 
as it is true for the classical case. 

Theorem 5.3 Assume that (r,/x) satisfies (po) (ind (VD). Then for any F e Wi 

NDLE (F) =^ LE (F) . 
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Proof. A modified version of Aronson's chaining argument gives the statement. The 
proof uses varying radii for the chain of balls. We give the idea of the modification (the 
other technical details can be seen following [34] or [14]). 

Denote 6 the constant in NDLE (F) and let 1 > 5' > arbitrary. If d (x, y) < Sf {x, n) 
the statement follows from NDLE, if 5'n < d{x,y) <n it follows from (po)- 

Assume that Sf {x, n) < d (x, y) < 5'n. Consider a shortest path tt between x and y, 
denote d — d{x,y) , 

R = f {x,n), S = f {y,n), A = B {x,d + R) U B {y,d + S). Let oi = x and 

ri = \Scof {oi,m)] 
and choose 02 G tt : d (oi, 02) = ri — 1 and recursively 

n = \Scof {oi,m)] (5.9) 

and Oj e TT : d{oi,Oi+i) = ri — 1 and (/(?/, Oj+i) < d{y,Oi). Denote Bi = B{oi,ri). The 
iteration ends for the first j for which y e Bj. Prom F e Wq and z^+i e B^ it follows that 

ci < 4r^^, < c, (5.10) 

and the from triangle inequality it is evident that 

d {zi, Zi+i) < 2ri + Tj+i < ( 2 + — ) Scof {zi, m) . (5.11) 



Here we specify Cq— {2-\- l/ci)~^. Let us recall the definition oil = 1 {n, d, A) 



taking the inverse one obtains: 



ft f d 

->maxCFL2, -), (5.12) 



Let us choose C in (5.12) (using F e l^i) such that 



By the definition of j 

d>J2^i>{j-l)j, 



d 
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consequently, j — 1 < I- 

[X,Zi) {Zi,Z2) ...Pm{Zj-i,y) . 



Now we use NDLE to obtain 



V / ^. V{xJ{x,m)) V{zj^i,f{zj^i,m)) 



^ ■ ■ j-i V(oi,ri) y (oj_i,rj_i) II (y) 

> mm . . . mm & — . . . — — 

~ z2eBi "' zj_ieBj_2 V{x, f {x, m))"'V{zj^2, f {zj-2,'m)) V {z^-i, f {zj^i,m)) 



> min ... min cP ^: 



z2eBi "' zj-ieBj_2 V {x, f {x, m)) V {z2, f {z2, m)) " V {zj-i, f {zj^i, m)) 
If we use (5.9) , (5.10) and {VD) it follows that 

c^'^Vd/) V{oi,ri) V{oj^2,rj-2) 



> min ... min 



> 
> 
> 



z2GBi z,_ies^_2 V {x, f {x, m)) y ^^^^ " V 

cxp [-C (j - 1)] 
exp [-CI] 



^.14) 



V 


{x, f {x, 


m)) 




cn{y) 




V 


f (yX, 


n)) 




cn{y) 





V{x,f{x,n)) 

From Lemma 13.6 of [14] we know that there is a c > such that 

I 



(5.15) 



Pn>C 



n—lm 



P 



if n > Zm + / — 1. Let us note that from (5.8) it follows that n — lm + l<'il which results in 



Pn {x, y) > c 



,n—lm 



^) {x,y) > 



J c^V (y) 



■ exp {-CI) 



> 



■•/' (y) 



V{x,f{x,n)) 



exp 



-C 



V(x,f{x,n)) 

F{:r.d i,r,y))\^ 



n 



This finishes the proof of the lower estimate. ■ 

Proof of Theorem 5.1. First of all we have seen that under conditions (VD) + (H) + (*) 
we have that 

E eWo or fw e Wo 

and {EE). If in addition G Ki, then both functions belong to Wi. On the other hand pv G 
Vi implies {aDpv) which is in the set of equivalent conditions (*), furthermore W^o^ = Wi 
which shows the equivalence of 2 and 3. Based on these observations the equivalence of 1 — 4 
and 5 is established by Theorem 4.1. The equivalence 5 and 7 is given in Theorem 5.2. 
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The implication 4 =^ 6 follows from Theorem 4.1 and 5.3. The reverse implication with 

respect to the upper estimate is also covered by Theorem 4.1 as well. For any F G Wi 
LE {F) implies (VD). This can be seen exactly as it is proved for F {x, R) = R^. The proof 
of UE (F) + LE (F) =^ PLE (F) can be reproduced following the steps of the proof of 
Lemma 8.3 in [34]. This shows the equivalence of 4 and 6 and proves the whole statement. 



Remark 5.3 Let us note that we can prove a slightly better upper and lower estimates (which 
are, in fact, equivalent to the ones presented). Denote d = d{x,y). Following the proof of 
the upper estimate in [36] (see the proof of Theorem 3.14 Remark 3.4) one can see that 

The intermediate estimate (5.15) gives a stronger lower bound: 

Pn{x,y)>— — — ^exp[-Cl{x,n,A)], (5.17) 

V{x,f{x,n)) 

where A = B {x,d {x, y) + f (x, n)) U B {y,d {x, y) + f {y, n)) ,n > d {x,y) . 

Remark 5.4 It is not immediate, but is elementary to deduce from (5.16) and (5.17) a 
spacial case of Theorem 5.1 if 

E (x, R)c^F {R) 

or 

p (x, i?, 2R) V (x, R, 2R) ^F{R). 

Such a result is presented in [36]. The key observation is that under (po) , (YD) , {H) the 
condition {E) implies /?'>!. The statements 1 — 5 and 7 of Theorem 5.1 are immediate, 
the two-sided heat kernel estimate 

exp [-Cm (n, d {x, y))] exp [-cm (n, d {x, y))] 
^ TTf — <Pn{x,y)<C — — — — (5.18) 

needs some preparation (here f (n) is the inverse of F {R) again). It follows from (5.16) and 
(5.17) and from the fact that for any fixed Cj > 0, x e F 

{Ciu, C2R) ^ Ix (Cgn, C4R) ~ m (Cgn, CqR) . 

In the very particular case when E{x,R) ~ R'^ one recovers from (5.18) the sub-Gaussian 
estimate (1.2) which is usual for the simplest fractal like graphs. 

List of the main conditions 
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shortcut 


equation 


name 


(Po) 


(1.5) 


controlled weights condition 


(VD) 


(1.3) 


volume doubling property 


{TC) 


(1.7) 


time comparison principle 


{ER) 


(3.1) 


the Einstein relation 


{MV) 


(2.16) 


mean value inequality 


DUE (E) 


(4.3) 


diagonal upper estimate 


DLE (F) 


(5.6) 


diagonal lower estimate 




(1.16) + (1.17) 


bounds on g 




(1.9) 


elliptic Harnack inequality 


UE {F) 


(1.18) 


upper estimate w.r.t. F 


PLE (F) 


(1.19) 


particular lower estimate 


TV T T — \ 7" 7 — I / T — 1\ 

NDLE [F) 


/a a\ 

(4.4) 


near diagonal lower estimate 


LE (F) 


(1.20) 


lower estimate 


PMV (F) 


(1.14) 


parabolic mean value inequality 


PSMV (F) 


(1.15) 


parabolic super mean value inequality 


(E) , (pv) 


(2.14) (2.4) 


E (x, R) or pv is uniform in x 


(aDpv) 


(2.7) 


anti-doubling for pv 


PH{F) 


(1.14) 


parabolic Harnack inequality 


RLEF 


(2.6) 


resistance lower estimate 
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